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Abstract 

We study non-perturbative real time correlation functions at finite temperature. In 
order to see whether the classical term gives a good approximation in the high temper- 
ature limit T ^ fiuj^ we consider the first quantum corrections. We find that for the 
simplest non-trivial case, the quantum mechanical anharmonic oscillator, the classical 
result is reliable only for moderately large times: after some time t* the classical ap- 
proximation breaks down even at high temperatures. Moreover, the result for the first 
quantum corrections cannot, in general, be reproduced by modifying the parameters 
of the classical theory. 
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1 Introduction 



Real time processes at finite temperature play an essential role in the physics of the early 
universe and of heavy ion collisions. A key quantity in scenarios of baryogenesis 
is the rate for electroweak baryon number violation (the sphaleron rate). In the broken 
phase the sphaleron rate can be computed with semiclassical methods |1|, ||, Q but in 
the symmetric phase ^ they are not reliable. Unfortunately, a direct non-perturbative 
lattice determination of the hot sphaleron rate is not available, either. 

The most promising approach to this problem ||^ is to compute the sphaleron rate 
in a classical real time simulation since the relevant thermal transitions are essentially 
classical. Considerable work has been done in this direction [[7|-[ll3|. 



Treating the dynamics of a classical gauge field system one is nevertheless faced 
with severe difficulties [|n|-[]l9|. The high momentum modes with k>T which do not 
behave classically, do not decouple from the dynamics. In general, these modes lead to 
ultraviolet divergences in the classical correlation functions which cannot be removed 
by introducing local counterterms in the classical theory [|1^, p!9[] . 

There is another question related to the classical approach which has hardly been 
considered so far: under which conditions is the classical approximation for the low 
momentum modes reliable? One systematic way of investigating the validity of the 
classical approximation is to compute the first quantum corrections in the ^-expansion. 
So far, the expressions have been derived only for quantum mechanics and scalar field 
theories pO[. However, these simple cases should already teach us something in spite 
of the fact that topological observables and the associated rate do not exist. In these 



models relevant observables might be related for instance to the damping rate [12, 13 



The purpose of the present paper is to evaluate the quantum corrections in the 
simplest non-trivial case, the quantum mechanical anharmonic oscillator. This study 
serves to estimate the feasibility of similar studies in field theories. Moreover, we believe 
that some of the general results might be carried over to that context. 

We find that while at small times the classical approximation is reliable, it breaks 
down at large enough times. The reason is that the functional form of the quantum 
corrections is qualitatively different from that of the classical answer, in a way which 
cannot be accounted for by modifying the parameters of the classical result. 

The paper is organized as follows. In Sec. ^ we discuss the formulation of the problem. 
In Sec. § we briefly discuss the harmonic oscillator and in Sec. ^ the anharmonic 
oscillator. The "symmetric" and "broken" cases of the latter are analyzed in more 
detail in Sees. ^, ^ and we conclude in Sec. ^ 
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2 The formulation of the problem 

We consider one bosonic degree of freedom q with conjugate momentum p and the 
Hamiltonian 

(2.1) 



where 



U{q) 



1, ,2^2 I 1 ^2^4 



(2.2) 



^ -2^; q + q + ^ 

We refer to the two cases of a positive and of a negative quadratic term as the symmetric 
and the broken case, respectively. Quantum mechanical (Heisenberg) operators are 
denoted by capital letters, for example 



Q{t) 



g(0)e 



h 



The finite temperature correlator we consider is 



C{t) 



1 r 



Q{t)Q{0) + Q{0)Q{t) 



ReTr 



-f3H(P,Q) 



Q(t)Q(0) 



relevant for the time dependence of 



Q(t)-g(o) 



(2.3) 



(2.4) 



(2.5) 



Here Z = Ti [exp{—i3H)] and (3 is the inverse temperature. Note that C{t) is an even 
function of t. 



In [^, the expansion 

C{t) = Cci(t) + Cn{t) + C,2(t) + 0{h^) 
was derived for C{t). The classical result is pT| , |22| , |15 



Ccl(t) 



dpdq 
2-Kh ^ 



(2.6) 



(2.7) 



where = J ^^e~^^^'P''^^ and qdt) is the solution of the classical equations of motion 
with the initial conditions qdO) = q, gc(0) = p. This expression corresponds to the 
prescription suggested by Grigoriev and Rubakov p. 

As for the quantum corrections, the contribution Cn{t) vanishes. The result to order 



is then 20 



C{t) 



-7 



dpdq 
2'nh ' 



-l3H{p,q) 



2/32 



24 



U"{q) + ^ (dl + U"{q)c)^ 



24 



g?c(it) 



-^g f^dt'U"'{qdt')){qdt'),qM},^ + 0{h^), (2.8) 
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where {, } denotes the Poisson bracket 



{f,9} = dpfdqQ - dpgdj, 



(2.9) 



and 

{f,9}o = 9, {/,^}n+l = {/,{/, 
Similarly, the expression for Z to order fi^ is 



(2.10) 



dpdq 



2/32 



24 



f2.111 



There are thus three kinds of terms in the h -correction to C{t), denoted by C^2(t), 
= a, 6, c: 



where 



24 



dpdq 
2'Kh ' 
_^ fTi^PX f dpdq 

dpdq 



-PiiMU"{q) 



Cci(t) - qqc{t) 



24 
24 



?gc(i), 



(2.12) 

(2.13) 
(2.14) 



-PHMq dt'U"'{q,{t')){q,{t'),q,{t)h. (2.15) 

ZTTtl Jo 



The termC^2^(t) is a sum of the correction to the partition function when it combines 
with the classical result Cci{t), and of the corresponding term in the numerator of 
eq. (^78|) . Eqs. ( p.l3| )- (|2.15| ) are the corrections we will evaluate below. 

One of the key issues of the present problem is the following: In the case of static 
time-independent correlators, it is possible (in a weakly coupled theory) to reproduce 
the results of the full quantum theory from a classical theory with a high accuracy, 
provided that the parameters of the classical theory are modified appropriately. This 
is called dimensional reduction [^, The question is then whether such a resum- 
mation might also work in the time-dependent case. Indeed, it has been proved that 
the resummation used in the time-independent context is sufficient for making the 
time-dependent two-point function in the scalar 0^ theory finite to two-loop order in 
perturbation theory and even for giving the corresponding damping rate the right lead- 
ing order numerical value . General arguments in the same direction were also given 
in |T^. The expansion in eq. ( |2.8D is, in contrast, non-perturbative: each term involves 
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contributions from all orders in the coupling constant. Let us therefore discuss the ef- 
fects of the resummation in the present context (see also [0). Of course, the problem 
of divergences does not occur unlike in field theory. 

First, consider dimensional reduction. Let us take as an example the "symmetric 
case" anharmonic oscillator, 

U{q) = \u:'q' + -/q\ (2.16) 

The starting point is then a 1-dimensional Euclidean field theory defined by 

111 r ^ rf^^ 

C = -{drqf + -u'q' + -/q\ Z = J Vqexp{-- drC). (2.17) 

According to dimensional reduction, this can be written as 

Z = const. X J dqoexp{—Scs), (2.18) 

where go is the zero Matsubara mode. The parameters in Ses are modified by the 
non-zero modes. The non-zero mode propagator is 

/ \ ^n+m,0 fry 1 n\ 

To order (which is a good approximation as long as (3hu) <^ vr), one can then easily 
calculate how the mass parameter in the effective theory is modified: 



, , 1 o 1 



^0 {27rnT/hy 4 



The change in the coupling constant is of order h'^ and thus does not contribute in the 
present ^^-calculation. 

Consider, on the other hand, eqs. (^), (|2ll| ). In eq. ( pTl| ), U" = + 3g^q^. The 
constant cu^-part of this expression does not contribute in eq. ( p.8| ) since it is cancelled 
by a similar part in the numerator, see eq. ( [2.131 ). The g^-part, on the other hand, 
can be exactly reproduced by calculating the classical partition function Z^i with cj^ 
modified according to eq. (p.20|): 



exp ( -pIu;V) ( 1 - K^3g\'] = exp (-(^lu'^.q' ] + 0{h'). (2.21) 



2 ^ J \ 24""/ "V 2 

Similarly, the — /3^[/"(g)-term in the square brackets in eq. ( p.SD is accounted for by 
the change in u"^ according to eq. (|2.2CI|) . Thus the term Cj^\t) in eq. (|2.13| ) is directly 
related to changing the parameters of the classical theory. However, there remain the 
terms C^2{t) and C]^}{t). On the other hand, qc{t) is still a solution to the original 
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Hamilton equations of motion. Hence the question is whether the ^^-effects can be 
taken into account by determining qdt) form the equations of motion with the modified 
parameter tu^g rather than u"^. This issue will be discussed below and we find that, in 
general, such a resummation does not take place. 

Finally, it should be noted that in the field theory case one is usually interested in 
a "rate" observable: a time independent constant determining the time dependence of 
some Green's function, for example the sphaleron rate or the damping rate. We are not 
aware of such an observable related to C{t) in the present context. We thus consider 
the general large-time functional behaviour of C{t). 



3 Harmonic oscillator 



In order to show in a simple setting how the ^-expansion works and to see what the 
structure of the perturbative solution is, let us start by considering briefly the harmonic 
oscillator. The classical Hamiltonian is 



H 



2^ 2 ^ 



In this trivial case, the correlation function in eq 
the result 

h 



Cit) 



(3huo 



■ ( tanh 
2uj\ 2 



-1 



(3.1) 

can be calculated exactly, with 

(3.2) 



Expanding in Ti, one gets 

Cel(t)+Q.(t) 



coscut 



12 



(3.3) 



The fact that it is the symmetric combination of Q{t)Q{fi) which appears in eq. ( p.4|) , 
removes the term linear in % from the result. It is seen that the quantum corrections 
change the amplitude of Cci(^), but not the frequency since u) is independent of energy. 
The classical ^"^-term is reliable in the limit jShuj <C 1, that is, at high temperatures. At 
low temperatures, in contrast, the T = result (with tanh = 1 in eq. ( p.2|) ) is reliable. 

How is this result reproduced by eqs. ( p.7|) , (|2.8|) ? The solution of the classical 
equations of motion is 

(3.4) 



Qcit) = q cosut H smutt. 



Substituting this into eq. (p.7[), one sees that the term proportional to p in qdt) does 
not contribute due to antisymmetry in p, and one gets directly 



Cci(t) 



1 cosut COSiUt 



(3.5) 
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where it was used that Zc\ = {Phuj)~^. 

As for the quantum corrections, the last term in eq. (2^) is proportional to the third 
derivative of the potential and thus does not contribute, C^2{t) = 0. The term C^\t) 
in eq. (|2.13|) does not contribute either, since U"{q) is just a constant. There remains 
a contribution from Cj^2{t) in eq. (|2.14|) , reproducing eq. 



4 Anharmonic oscillator 

Let us then move to the less trivial case of the anharmonic oscillator. Here and in the 
following we use u, g and 

y«-$ (4.1) 

to introduce the dimensionless variables q,p,t, (3, E: 

q = -q. P=^P, i = iut, P = /3Vo, E = ^. (4.2) 

This rescaling serves to show the parameter dependence of the final non-perturbative 
result more clearly. At the same time, it makes the coupling constant equal to unity 
so that if one wants to compare with perturbation theory, one should go back to the 



original variables. In terms of the rescaled variables the potential in eq. ( |2.2|) reads 



\ _ j Vo(g'^ + 2g^), the "symmetric" case (a'W 
K)(g'-1)', the "broken" case ' ^ ' 

A dimensionless combination to which 7i can be attached is 

-4^. (4.4) 

The quantity naively governing the semiclassical expansion is hence e^. This is multi- 
plied by some dimensionless function /(/3, i) which may scale approximately with some 
power of /3 for given t. For instance, in the case of the harmonic oscillator, f{P,i) 
scales as so that the real expansion parameter is 

(e/3)2 oc {f3hujf. (4.5) 

One of the issues below is how the function f0,i) behaves in the anharmonic case as 
a function of (3. 

With the rescaling performed, one can also write C{t) in a dimensionless form. 
Factoring out the scale u'^/g'^, the classical correlation function is 
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where 



^cl 



Here we utilized the symmetry of the integrand in g ^ —(j,P 
the quantum corrections in a dimensionless form, 

,2\ 



dp I dqe-^^. (4.7) 
—p. We can also write 



UJ 



(4i 



We will then discuss the "symmetric" and "broken" cases separately. 



5 The symmetric case 
5.1 Numerical results 

The detailed form of the classical solution qdt) and of the integrals appearing in the 
symmetric case is discussed in Appendix A| . The expressions to be evaluated are in 



eqs. (|4.6|) , ( |A.9|) -( [An^ ). We have done the evaluation numerically, as well as analyt- 
ically in certain regimes of Let us discuss the numerical result first. The curves 
are displayed in Figs. 

The qualitative features of the solution are the following: Both the classical solution 
Cci(t) and the quantum correction C^\t) + C^2{t) approach zero at large times. The 
time scale it takes for the amplitude to diminish depends on /3, being roughly propor- 
tional to l3, and being somewhat larger for the quantum corrections. The reason for the 
attenuation is the destructive interference of the continuum of classical solutions with 
different frequencies. This feature does not fully persist in the quantum case where the 
energy levels are discrete: rather the behaviour is "almost periodic" I^T] on a larger 
time scale. Indeed, already the term C^2 (t) in Fig. ^ behaves in a manner qualitatively 
different from C^\t), cj^2{t): it has a constant amplitude at large times. 

Let us discuss these features and their implications in more quantitative terms. 



5.2 The large time limit 

We are mainly interested in the large time behaviour of the correlation function. Ordi- 
nary perturbation theory breaks down for large times and is therefore excluded. This is 
due to the secular terms in the perturbative series: at lowest order the solution to the 
classical equations of motion is proportional to cos^ut+a) while the next order contains 
a term proportional to g'^tsm{ujt + a). Thus, by dimensional analysis, straightforward 
perturbation theory only works for 

i = ujt<^—(3 = Wo- (5.1) 
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Figure 1: The classical correlator Cci(t) in the symmetric case. The thin line represents 
the analytic approximation of Sec. [5.2| . The difference between the thin and thick 
lines in the regime 1 < cut < /3Vo is due to higher order perturbative corrections in 

The way to avoid the secular terms is to use the exact frequency Vt{E) inside the 
trigonometric functions appearing in the perturbative series. The perturbative series 
for the classical solution of eq. ( |A.3| ) is obtained from eq. (|C.6|) . In the phase space 
integration one then has to compute (after a change of variables according to eq. ( |A.15| )) 
the dimensionless integrals 

Jnm{P,t) = (^) / dEe-''''E^sm{mQ{E)t). (5.2) 

JO 

In general, this is difficult to do. Fortunately, an exact evaluation is not necessary if 
one is interested in the large time limit ut ^ 1. Then it is sufficient to keep only the 
first two terms of the low energy expansion of the exact frequency, 

Q{E) =uj[l + ^ciE + 0{E^)], (5.3) 
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Figure 2: The quantum correction C^\t) + C'^2{t) in the symmetric case. We have 



divided out the naive expansion parameter (e/3Vo) 



where E = E/Vq and Ci = 3/4. In this approximation we find 

cos(a;t + (n + l)ip) 



JnlW,t) 



{n+l)/2 ' 



sin(ci;t + (ri + 1)^9) 



where 



arcsm 




CiUt 



(5.4) 



(5.5) 



From these expressions one can see that in the region ut^ ^(3 the terms which have 



been neglected in eq. ( |5.4| ) are suppressed by at least one power of l/{ujt). The reason 
is that each power of E in the phase space integrand of eq. (|5.2| ) gives one power of 
ll{ut). In the region 1 <^ ut <^ [3Vq, the terms neglected are suppressed by g^fijuj'^^ 
corresponding to higher order perturbative corrections. Note that the approximation 
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Figure 3: The quantum correction C^Jo(t) in the symmetric case, compared with the 



analytic approximation of Sec. 5.2 



( ^.4|) is vahd also for uot ^ /3Vo where the perturbative expansion breaks down. Thus 
the large time expansion for C{t) can be obtained from the low energy expansion of 
the phase space integrand. 

Using this expansion, we find for Cc\{t) for ujt^ 1, 

,5 



UJ 



(5.6) 



If /5 = ^ 1, there is an overlap of the "perturbative region" of eq. ( |5.1| ) and the 



large time region: for moderately large times 1 ^ cut -C we recover the leading 
order perturbative result. If cut ^ in contrast, eq. (|5.6| ) simplifies to 



16 1 id cos LUt 



(5.7) 



9 nZci (ujty 

That is, for large times, the classical correlation function oscillates with the 'tree level' 
frequency u and with an amplitude which decreases as 1/t^. Comparing eq. ( p.7| ) with 
the corresponding result for the harmonic oscillator, eq. (|3.5|), we see that eq. 



IS 
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non-perturbative since its functional form cannot be obtained by adding corrections 
multiplied by positive powers of to the harmonic oscillator result. 

Let us now note that if a resummation according to eq. ( |2.20|) would take place, then 
the 'h? quantum result should be obtained by replacing uj"^ ujI^ in eq. ( ^.7|) , that is 

^resummed , . ^^^^/^l ^ , 2 1 U'^ (3^^^ Sinut 



Cci(t) + -, (5.8) 

9 nZc\ g'^ ojt 



where hi is some number. We show below that the true Cf^2{t) is not of the form in 
eq. (§3). 

We start with C^^(t). It was pointed out already in Sec. Q that this term is related 
to the replacement u"^ — > cj^^ in the Hamiltonian appearing in the Boltzmann factor. 
To be more specific, the term uj"^ in U"{q) cancels in eq. ( ^.1^ ) and in the limit tut ^ 1 
we find 

d^\t)^lihg(3nq\,C^,it). (5.9) 

The contribution proportional to {q^qc{t))c\, on the other hand, has one additional 
power of E in the phase space integrand compared with the classical case and is thus 
suppressed by a factor l/{ut). From eq. (|5l9|) it is obvious that the quantum correction 



C^\t) shows the qualitative behaviour indicated in the first term in eq. (|578|) : it is 
small compared with the classical result if (3huj <^ 1 and this holds even for arbitrarily 
large times. 

Next we consider the quantum corrections containing the derivatives d"^, which 
we have denoted by C^^}{t). These derivatives acting on the trigonometric functions in 
qqdt) give extra factors of t. When expanding the integrand in powers of energy one 
has to count t as E~^. For cut 1 we find 

C^Ht) ^ {^MM - 9u;tMf3,t) - ^-{ujtfM(3,t)^ . (5.10) 

The individual terms in the curly brackets behave as sin(c(jt) / (ut) for large times, which 
is the expected behaviour in eq. (|5l8| ). Such a result would at the same time indicate 
that without resummation, the semiclassical expansion breaks down for 



when the correction term in eq. ( |5.8D is as large as the leading term. However, we find 
that this does not occur: in the limit tut ^ 1 the individual terms in the curly brackets 
in eq. (|5.10|) cancel at leading order in l/(co't). Therefore the amplitude of C^^it) 



decreases as l/(u;t)^ for large times. Thus C^2^(t) is small compared with the classical 
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result at high temperatures. The corresponding suppression factor, however, is not 
in general given by {f3huj)'^. There are terms proportional to having different 

dependences on the temperature: expanding eq. (|5.1CI|) gives terms oc (3^ while the 
subleading terms in the low energy expansion are proportional to (3. We have not 
calculated these terms analytically. The numerical result for the sum of Cll^\t) and 
C^2{t) is shown in Fig. ||. 

Finally we consider the correction cj^2{t). The result for ut ^ 1 is 

^n^^t) ^ ^j;^-M)"{Mf^^t) - jujtjn{f3,t)}, (5.12) 

,4 



9" 

1 1 ri^ 



which for cut ^ becomes 

~ - — coswt. (5.13) 

Thus at large times C^^}{t) oscillates with the "tree level frequency" but with a time 
independent amplitude. This behaviour is qualitatively different from the classical 



case. Comparing eqs. (|5.7|), ( |5.13| ) we see that Q2^(t) becomes as large as the classical 
correlator for t ~ t^, where 

UJ^ 1 

uU = ^ = -, (5.14) 

and for t > t^, the semiclassical approximation breaks down. 

The correction in eq. ( |5.13|) is clearly not of the form allowed by eq. (|5.8|) . Since 



there is a term of a functional form not allowed and the allowed sin u;t/(co't)-term does 
not emerge, we conclude that a resummation according to eq. ( p.20| ) does not take 
place in the large time limit. Neither can one understand the result as a resummation 



with a correction factor different from that in eq. ( 2.20 ). Since a resummation cannot 
be made, the semiclassical expansion breaks down at the time given by eq. (|5.14| ). 

It can be checked from Fig. ^ that for cot > 10 the analytic approximation for C^^}{t) 
indeed gives quite an accurate estimate of the exact numerical result. 

To conclude, let us point out that the qualitative features found, together with the 



"almost periodic" behaviour at time scales t>t*, can be reproduced with the 
following approximation. Writing the full quantum result in eq. ( p.4|) in the energy 
basis, one gets 

C{t) = Z-iRe^e-^^™e^*(^™-^")|(m|Q|n)|2. (5.15) 



m,n 



Approximating the energy levels to first order in g 



2 



E^ = hco[n+^^ + ^^(n' + n + ^], (5.16) 
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and the eigenstates to zeroth order, one gets 

Tn=0 



C{t) ~ Z-'^ E (m + 1) (e-^^-+^ + e-'^^™) cosf ^'"+^^ t. (5.17) 



The behaviour of this solution for small e = Tig^juj^ follows the classical solution in 
Fig. |I] until the time scale is of order t ~ = 4/e, but then the periodicity sets in so 
that at the time scale t ^ 8t^,, the structure around t = in the classical solution is 
repeated. This is the reason for the breakdown of the classical approximation. 

6 The broken case 
6.1 Preliminaries 

In the broken case, the classical Hamiltonian is 



H = Vn 



2p^ + if - if . (6.1) 



There exists, of course, an enormous literature on this system. In the present finite 
temperature context, it has been previously studied by Dolan and Kiskis and by 



Bochkarev 22 



Quite a lot is known about the qualitative behaviour of C{t). In general, the solution 
can be written as in eq. ( |5.15|) . Since the solution is a sum of periodic contributions 



corresponding to the different energy levels that can be excited, C{t) is "almost peri- 
odic" In particular, the lowest frequency appearing is determined by 



3 



AE = Ei-Eo(x huexp 1 , (6.2) 



implying that the symmetry is restored already at T = in the sense that the 
correlator averaged over a long enough time period vanishes. In contrast, the classical 
result Cci(t) has a non-zero limiting value for t ^ oo, in which the symmetry is only 



partially restored and all the oscillations die out ^T|. The oscillations die out, like in 
the symmetric case, due to the destructive interference of the continuum of classical 
solutions with different frequencies. 

The fact alone that the classical result does not show the expected qualitative be- 
haviour of the full result, indicates that the classical result is not generically applicable. 
We study this problem in more concrete terms below by evaluating the ^^-corrections. 

Note that, in seeming contrast to what was just pointed out, the system in eq. (Ol) 



has also been used to illustrate that the classical approximation is applicable to some 
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real time problems 0. The reason for the difference is that the situation we con- 
sider is different from the one in [Q: we have a strict equilibrium situation at a finite 
temperature which is also what is considered in ||21], and which occurs in the 
real time sphaleron rate simulations. The consideration in 0, in contrast, concerns 
a non-equilibrium symmetry-restoring rate obtained by taking an initial state where 
the system is prepared in one of the minima. In the strict equilibrium case, one can- 
not define such a rate. Still, the problem of the general applicability of the classical 
approximation to real-time problems remains. 



6.2 Numerical Results 



The form of the classical solution qdt) for the broken case is discussed in [Appendix B 



The numerically evaluated classical correlator Cc\(t) is shown in Fig. ||, and the quan- 
tum correction C^'' {t) + C^2 {t) in Fig. |. 

The most notable difference with respect to the symmetric case is that there is 
a constant part in the broken case results. The energy integrand for Cci(t) is for 
illustration shown in Fig. ^ where the emergence of the constant part from E < 1 
can be seen. It is evident from Fig. ^ that the partial symmetry restoration in the 
classical result is the stronger the higher the temperature is and from Fig. |^ that 



there is a further symmetry restoring effect from the quantum corrections. It is seen 
in Fig. ^ that at high temperatures (/5Vo ~ 0.5 — 2.0) the quantum corrections are 
roughly proportional to the naive expansion parameter (e/5Vo)^ = {jPhu)'^ which has 
been factored out. 

Let us then discuss C^2 0,i). Its numerical evaluation turns out to be very difficult 
for large i. The reason is that the energy-integrand is highly peaked and oscillatory 
around unity. To see this, note first that at t = 0, the integrand in eq. ( |A.12| ) vanishes. 
Moreover, the integrand involves terms ~ sm^l{E)i, in analogy with eq. ( |6.8| ) below. 
Hence a particular energy region will contribute provided that 

^iE)i>l. (6.3) 



Let y = \E - 1\ <^ 1. Since K{k) ~ ln(4/A;') close to k = 1 {k' = Vl - A;^), one gets 
from eq. ( [B.4| ) that 

{ 27r E ^ 1. 

ln{64/s,)' _ ^g4) 

ln(64/s/) ' ~ 

Eq. ( |6.3| ) shows then that the energy-integrand can be large in the region 

y>e-'. (6.5) 
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Figure 4: The classical correlator Cci(t) in the broken case (thick lines), together with 
the analytic approximation of Sec. B]3 (thin lines). 



On the other hand, the second partial derivative in eq. ( [A. 121 ) will involve 

dln{E) = {dpEfdln{E) + . . . , (6.6) 
where dpE = 4p. Hence according to eqs. (p. 41), (|6.5D, 



,2 



~-<e'*- (6.7) 

y 



Thus the height of the peaks around E = 1 grows exponentially with time, and the 
peaks move closer to E = 1. The width of the peaks is diminishing, but their height is 
growing faster so that they give an increasing contribution. In fact, the highest peak's 
contribution from ^ < 1 (where the peak gives a positive contribution) and from 
E > 1 (where it gives a negative one) to a large extent cancel, but the cancellation is 
not complete and one has to account for it very precisely in the numerics to get the 
remaining contribution correctly. This is why we cannot go to large i. In practice, we 
can reliably calculate C*^2'*(/3, i) only up to t = 15, when the highest peaks in the energy 
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Figure 5: The quantum correction C^^ (t) + (t) in the broken case. 

integrand are of height ~ 10^° (for (3Vo ~ 2) at 6E ~ 10^^ around unity. The energy 
integrand is shown in Fig. ^ and the result of the integration in Fig. |^. 

6.3 The large time limit 

Consider first the classical correlation function. The form of the solutions in eq. ( |B.2| ) 
can be read off from eq. ( |C6|) . It is seen that for E < 1, qdt) contains a constant 
part in addition to the cosines. The 0-integral obtained with the change of variables 
in eq. ( [A.15|) , gives then 



rK[k) ^ 2 |- oo 2n n 

The cosines in eq. (|6.8|) give contributions which vanish in the limit t ^ oo as in 
eq. (|5.4| ), see below. Hence one gets from eqs. ( |4.6|) , (|A. 15|) that the constant part 
surviving is 
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Figure 6: The energy integrands for Cci(^) and C^2{t) in the broken symmetry case 
for l3Vo = 2. In the hmit — > oo only the region E/Vo < 1 contributes in Cci{t). 
The energy integrand for C^^it) has been shown on a logarithmic scale. It involves 
essentially the second derivative of the classical integrand, which is why it has very 
high peaks (~ 10^° already at ut ~ 15) around E/Vq — 1. 



One may also try to compute the time dependent part for ujt ^ 1 in the same way 
as in the symmetric case. There we saw that the limiting behaviour for large times can 
be obtained from a suitable low energy expansion of the solution to the equations of 
motion. In the present case it is obvious that this expansion cannot be convergent when 
E approaches unity. One may argue, however, that for large times only the solutions 
with small energies are relevant and that this expansion still works. We find that 

C,i{t) = C,i{t ^ oo) + -L^^^^'d£;e-'^^(|icos[l](E)t] + 0{E')y (6.10) 

Replacing the upper integration limit by oo and keeping only the first two terms of the 
low energy expansion of 

fl{E) ^ V2u;{l - diE + 0{E^)) (6.11) 
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Figure 7: The quantum correction C^2{t) in the broken case, divided by (e/5Vo)^ = 

where di = 3/16, we obtain for cjt ^ 1, 

^ / X ^ / X 1 1 cos( V2cjt - 2(p) 

Cci t ^ C^iit - oo -f — — -— - J[ 6.12 



where now 



ip = arcsm = . (6. Id) 

,^{pW+{d^v^j 



It can be seen in Fig. ^ that eq. ( |6.12| ) is indeed a good approximation at large times. 

The integrals appearing in the quantum corrections C^\t) and C^2{t) are qualita- 
tively quite similar to that appearing in Cci(t). In particular, there is a constant part 
in these corrections which can be evaluated in the same way as eq. (|6.9| ) . It is seen that 
the constant part tends to further restore the symmetry compared with the classical 
result, see Fig. |[ 
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For the quantum correction C^2^(t), in contrast, the "small energy expansion" does 
not seem to be applicable. We have computed the solution but it does not agree with 
Fig. 1^. However, this need not be a surprise since, as discussed, it is not guaranteed 
that the small energy expansion works in the broken case due to the singular nature of 
the point E = 1: the energy integration extends beyond the radius of convergence of 
the small energy expansion. Moreover, the integrand in C^2 {t) is qualitatively different 
from that in Cci(t). A simple example where the small energy expansion would not 
work is given by 

m ^ f dEe~>^?}l^!^. (6.14, 

•''0 7l{l-E) 

At £ ^ oo the integrand makes a delta-function, and f{i oo) exp(— Yet an 
expansion in E of the denominator around E = and an integration term by term, 
gives a result which oscillates around zero. 

We could not find any other analytic way of evaluating the energy integral for C^2{t), 
either. The integrand is very complicated around E ^ 1. Thus we can only mention 
some general features of the solution. 

First, note that the numerical result in Fig. |^ shows that there is a growing negative 
contribution at large i in C^2{t)- This seems to arise from E a bit larger than unity. 
To estimate very roughly when this kind of a contribution can be important, note that 
then the peak heights must be such that the exponential suppression cannot hide them 
any more, that is 

e-V*>l. (6.15) 

Hence one starts to get an effect at t > /3. 

As to the functional form of the solution, it looks roughly like —i^ at large i. It is 
easy to see that a linear in t behaviour cannot occur, since it follows directly from the 
definition in eq. (|2.4| ) that C{t) is symmetric in t. For j3Vo = 0.5, 1.0 in which case the 
asymptotic behaviour is obtained earliest, the leading term of C^2{t) can be fitted at 
ci;t ~ 8 . . . 15 for instance as 



d'^it) ~ ^4^' -OM{ujt)^ . (6.16) 



The conclusions one can draw from the broken case seem rather similar to those from 
the symmetric case. The quantum correction C^a"* (t) behaves in a manner qualitatively 
different form what was observed for Cci(t). Moreover, the difference is such that it 
cannot be accounted for by a simple resummation of the mass parameter cj^. As the 
classical result in Fig. § is of order unity and the fit in eq. (|6.16| ) would suggest the 



behaviour e'^{ujt)'^ for the quantum correction, one would expect that the semiclassical 
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expansion breaks down at 

a;t~.t. = ^=(-) . (6.17) 

In eq. ( |5.14|) in the symmetric case it was rather observed that ut,, = 1/e. However, 
the fit in eq. ( p.l6| ) should not be taken very seriously as the interval is very small, 
and the main point is that the time scale for the breakdown seems to be determined 
by 1/e. 

Finally, let us point out that from the general form of eq. ( |5.15|) , one might have 
expected that at finite temperature the asymptotic values of C{t) are oscillating be- 
tween positive and negative values. At zero temperature the time scale would be 
~ exp[2v^/(3e)] according to eq. ( p.2|) . Thus the quantum correction C^2{t) seems 
to restore some of the qualitative features missing in Cci(t), in the sense that the be- 
haviour in Fig. ^ looks like the beginning of an oscillation with a large time scale. The 
difference from the zero temperature case, however, is that the time scale associated 
with the oscillation is not exponential. 



7 Summary and Conclusions 

We have studied the classical finite temperature real time two-point correlation function 
and its first quantum corrections for the anharmonic oscillator. The expansion around 
the classical limit is made in powers of h, so that each order contains all orders in the 
coupling constant g^. 

One can identify three different time scales in the results. In the symmetric case 
(Section |), these are 

~ 1, ujt P = —^(3, ut ut^ = - = — -. (7.1) 

As long as ut -C f3, perturbation theory works and the correlation function oscillates 
with period cut ~ 1. In the non-perturbative region ut > (3, the correlation function ap- 
proaches its asymptotic form. We have developed a large time expansion which allows 
to address also the time scales ut ^ (3. In this regime the amphtude of the oscillations 
in the classical result attenuates due to the destructive interference of solutions to the 
equations of motion with different energies. This attenuation cannot be associated with 
a damping rate. Finally, the time scale is associated with the quantum corrections 
and becomes infinity in the formal limit ^ 0. There is a hierarchy utt, ^ (3 provided 
that (3huj < 1. 
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The general result of our study is that at the non-perturbative time scales out > (3, the 
form of the quantum corrections differs qualitatively from that of the classical result. 
The semiclassical expansion breaks down at t ~ when the quantum corrections 
become as large as the classical result. Moreover, we found that these large corrections 
cannot be resummed by modifying the parameters of the classical theory. 

On the other hand, the first quantum corrections to the classical correlation function 
are small for i <^ t*. Prom this we would expect that in this region the classical limit 
gives a good approximation for the full quantum mechanical correlation function. The 
expansion parameter for the quantum corrections in this region is not just the naive 
one {PhuY, but e{(3hw) and appear, as well. 

An essential question is then which of the discussed features might be carried over to 
field theory. Unfortunately, we cannot say very much about this. However, certainly 
the present study does not encourage one to believe in the generic applicability of the 
classical approximation in the high temperature limit for time-dependent quantities 
at arbitrarily large times. On the other hand, there arc also obvious features which 
cannot hold in a four-dimcnsional field theory: for instance, we found that the time 

docs not depend on the temperature. This is unlikely to be true in the pure SU(2) 
theory, say; dimensionally, the classical time scale not involving h is {g'^T)~^ in that 
case and the time scale proportional to h^^ is {hg^T)~^ . 

It would be interesting to extend the present type of an analysis to field theory to 
be able to make more concrete conclusions. Unfortunately, a straightforward evalua- 
tion of the quantum correction C^2^(t) was numerically quite demanding even in the 
present case, in particular for the "broken" case where the modes with E /Vq ~ 1 are 
rather singular. In the field theory case, the partial derivatives of the classical solution 
with respect to the initial conditions would be replaced by functional derivatives, mak- 
ing things even more complicated. Still, one might hope that the scalar field theory 
analogue of the symmetric case would allow a non-perturbative investigation of the 
quantum corrections in the damping rate. 

Finally, let us point out that as it appears that the classical approximation does 
not describe the large time behaviour at least in the present case, it would perhaps 
be useful to consider the feasibility of other approaches. In principle the problem 
can be solved non-perturbatively using Euclidean simulations and spectral function 
techniques. The anharmonic oscillator considered in this paper might be a suitable toy 
model for developing techniques for such studies, since it appears that there is some 
non-trivial structure even in this case. 
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Appendix A 

In this appendix, we give some details concerning the classical solution qdt) and the 
computation of the quantum corrections to Cci(t) in the symmetric case. We use the 
rescaled variables defined in eq. 
The classical Hamiltonian is 



H = Vo 

Let us introduce some useful notation: 



2p' + f + 2q'. (A.l) 



+ ^ G(2,oo), k'^ = l-e<0, ~k = k-'<l, 



1+E-l 



qI^ = Vl + ^-1, ^ = ^. (A.2) 
The solution of the classical equations of motion is then of the form 

qdi) = gmaxCn^(A;wt + k(f)), (A.3) 



where cnk{u) is a Jacobi elliptic function (see [Appendix C| ). Here the initial conditions 
of qc{t) at t = are given by 

gc(0) = q = gmaxcn^(/i;0), q d^) =p = -gmax^A; sn^(A;0)dn^(A;0). (A.4) 

For g > these can be inverted to give 



sign(p)fcF( arcsin. ± — — — , 

^ \ 9'max 



l--^,k], (A.5) 



where F{(f), k) is the normal elliptic integral of the first kind (see [Appendix Q ). If g < 0, 
one should replace the result in eq. ( [A.5[ ) by 

^ 2sign(0)Reir(A;) -0, (A.6) 

where 
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(the case A; < 1 is relevant in |Appendix B| ) . Note also that the frequency of the classical 
solution depends now, in contrast to the harmonic case, on the energy: according to 
Appendix Q , the period of cn^ is 4:K{k) so that it follows from eq. (|A.3| ) that 



2tt 



71W 



P{E) 2kK{k) 



{A.i 



where P{E) is the period of qdi)- 

Given the classical solution, the classical correlation function Cci(t) can be calculated 
from eq. ( [4 .61) . The classical partition function of eq. ( [4.7|) is given by 



1 f¥ 



(A.9) 



where Ki,{x) is the modified Bessel function of the second kind. The quantum correc- 
tions in eq. (|4.8| ) are given by 



dp rfge-^V[4T'C'ci(t)-ggc(t)], 

JO 

POO « . 

dp dqe-^^\-{3f + 1) 



+2(3q'{q' + ly + 2f3p'{3q' + 1) qq,{t). 



(A.IO) 



(A.ll) 



To get eq. ( [A.ll|) partial integrations with respect to p and q were performed. 

Concerning C^2 0,i), it is useful to make a canonical transformation to variables at 
time t' to evaluate the Poisson bracket in eq. (|2.15| ), to change then the time integration 
variable and to perform one partial integration with respect to p. The result can be 
written as 



cl'':)0,i) = ^ r dp r dqe-^^ I ^ di' A^mS' - 1) - 



1 



(A.12) 



The partial derivatives of qdi') can be evaluated numerically, or even analytically using 
the formulas in [^. As an example, the first derivative is 



dp 



P 

2P 



2 k'^ ^ . 



+snrfdnrv 



wq 



E-,{v) - E-,{k<j>) wt 

k 



kk'"^ p 



^max 



wkqi,^^ q 
2 p. 



+ 1 



v=kwt+k(j> 



(A.13) 



where Ef^{v) is defined in eq. 
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Finally, note that since the variable appears in a simple manner in qdt) in eq. ( |A.3|) 



it is convenient to make a change of integration variables. We go first into the canonical 
action-angle variables (/,«), and then from these into energy E and the variable (p, 
using 

The integration measure can then be written as 

/oo I'OO POO „ rReK{k) 1 

dp dq[...]= dE d<p — [...]. (A.15) 

-oo Jo Jo J^ReK{k) Aw^ ^ 

Appendix B 

In this appendix, we describe the classical solution qdt) and formulas for the quantum 



corrections to Cci(t) in the broken case. The classical Hamiltonian is in eq. ( |6.1| ). 
In accordance with eq. ( |A.2| ), let us introduce some notation: 



= ^^€(0,2), k'^ = l-k\ k = k-\ 



1 + \JE 

?max = 1 + W=^^. (B.l) 



Then the classical solution is 



^ (t\ _ / ±'?maxdnfc(wt + 0), E <l 



where the relation to p, q is (for g > 0) 



-signfplF arcsinJ^ 1 - J— , E<1 



sign(p)A;F( arcsinWl — — ,k], E>1 



(B.3) 



Note that E < 1 corresponds to modes which do not cross the barrier, whereas the 
^ > 1 modes do cross it. The frequency of the classical solution again depends on 
energy and is given by 

27r f TTTTT, E < 1 



The classical partition function of eq. ( [4.7|) is 

Z.m = M r dqe~^(^'-'^' = \ [le-^/'[V27rh0/2) + K,0/2)], (B.5) 
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where /jy, K^y are modified Bessel functions. Tlie quantum corrections to tlie classical 
result of eq. (|4.6|) are 



oo ^oo 



C^\p,i) = I dp I dqe-^^f[Z^,'C,,{t) - qUi)], (B.6) 



-oo JO 

/■oo f _ _ 

32 / / j-„-l3E 



-P' / dp / rfge-^^ (1 - 

3 J-oo JO L 



+2/3q'{q' - 1)^ + 2/3p^(3g^ - 1) qUt) 



(B.7) 



The correction C^^}{i3,t) is given by eq. ( |A.12| ). The partial derivatives of qc{i') with 
respect to p can again be evaluated analytically. For example, for £' < 1 we get 



dp 



^(2 - k'^ fqms.^{ dufcM 



^2 2 
Q'max Sn^M 



+snfcMcnfc'U 



Eu{u) - Eu{(t)) 



1 P -^2 



+ 1 



k''^w q 



qp 



fB.81 



M=«)t-|-(/) 



and for E > 1 the expression is the one in eq. ( A.13| ) with the replacement gj 
-^max) except for the q^^^ appearing in gVCax- 

Finally, a change of integration variables can be made according to eq. (|A.15|) 



2 

max 



Appendix C 

We discuss here briefly some of the basic definitions of the Jacobi elliptic functions 
used. The notation follows P7[ . 

Let F{(f), k) be the normal elliptic integral of the first kind, 

da 



F{<P,k) 



'o vT— Psiii^ 

Then the complete elliptic integral of the first kind K[k) is defined by 



K{k)=F{l,k)= 



The associated nome is 



where k'"^ = 1 — k"^. Similarly, 



q = exp 



-IT- 



{1 - t^){l - kH^) 

K{k') 
K{k) 



famj; u 



da 



rv I f 

E{(l),k) = / (iav 1 — A;2 sin^ a, u= . _. 

JO JO vl — A;2 "=1" 



sm a 



EJu) 



E{a.mkU, k) 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



25 



define the normal elliptic integral of the second kind E{(j), k), the amplitude function 
am^M, and the function Ek{u). 

The Jacobi elliptic functions are denoted by dn^w, sn^u, cufcU. They are defined by 
(0 < M < K{k)) 

/•I dt _ dt 

Ln,u ^(l-t2)(^2_^/2) icn,« ^(1 - (^/2 + p^2) 

(C.5) 



^(l-t2)(l-A;2t2)' 



Note that sn^w = sin(amfc-u). The periodicity of dn^w is 2K{k), and that of sn^M, cn^-u 
is 4:K{k)] dufcM, cufcM are symmetric while sn^u is anti-symmetric. The functions ap- 
pearing in the classical solution qdt) for the anharmonic oscillator are dn^w and cn^w. 
We need their Fourier expansions, 



vr 



2K{k) 



q nnu 
cos 



dufc-u 

CUfcM 

where q is the nome in eq. 
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